Baseline amplicon sequencing data for the ITS2 region in the green honey of Banggi Island, Sabah.

Green honey, was discovered on Banggi Island, Sabah, showing high in essential amino acids and chlorophyll derivatives. Despite its lucrative market potential owing to its distinctive color, uncertainties persist regarding its nature. This study leverages amplicon sequencing by targeting micro- and macro-organisms present in honey environmental DNA (eDNA) using Internal Transcribed Spacer 2 (ITS2) region, enabling the identification of floral and microorganism sources that represent the honey's composition. The investigation into green honey from Banggi Island concerns the prevalence of honey adulteration and authenticity for economic gain. Adulteration methods, such as the addition of sugar syrups, compromise honey purity. Using a sequencing approach would help in determining the geographic origin and verifying the authenticity of the honey. The study aims to identify plant species or microorganisms in honey's eDNA. To authenticate honey, we utilized ITS2 with Illumina sequencing, exploring the diversity of green honey samples. Raw sequence reads obtained for the green honey sample revealed 1,438,627 raw reads, with a GC average of 49.22 %. A total of 44 amplicon sequence variances (ASVs) were identified, including three genera: Zygosaccharomyces with two species, Fraxinus with three species, and the genus Ficaria with only one species. Their respective relative abundances were 98.55%, 0.94%, and 0.51%. Zygosaccharomyces rouxii and Zygosaccharomyces mellis were identified as the pre-dominant yeast species in honey, while the Fraxinus and Ficaria genus represent common plant species in Sabah, particularly in Banggi Island. The dominance of Zygosaccharomyces species aligns with their known prevalence in honey, affirming the reliability of our findings. The presence of Fraxinus and Ficaria in the honey sample correlates with its abundance in the local environment. This amplicon sequencing approach not only contributes to our understanding of green honey composition but also serves as a valuable resource for authenticating honey origin in Malaysia, particularly for green honey from Banggi Island, Sabah. Our study pioneers the application of ITS2 amplicon sequencing for green honey amplicon sequencing, providing valuable insights into its composition and origin. This methodology, with a focus on eDNA, contributes to the authentication and quality determination of honey in Malaysia, addressing the pressing concerns of adulteration and variability in production practices.

Green honey of Banggi Island: A preliminary anti-diabetic study on zebrafish model.

Zebrafish is a developing vertebrate model with several advantages, including its small size, and high experimental efficiency. Malaysia exhibit one of the highest diabetes rates in the Western Pacific and incurring an annual cost of 600 million US dollars. The objective of the study is to determine the antidiabetic properties of green honey (GH) using a zebrafish model. Adult zebrafish, aged 3-4 months, were subjected to overfeeding and treated with streptozotocin (STZ) through intraperitoneal injection (IP) on days 7 and 9. The study assessed the oral sucrose tolerance test (OSTT) and the anti-diabetic effects of green honey. The evaluation was conducted at three time points: 30, 60, and 120 min after treatment and sucrose administration. The study utilised a model with a sample size of 5. The study was performed in six groups. These groups are (1) Normal control (non-diabetic, no intervention), (2) Normal control + GH (non-diabetic, supplemented with GH 3 µl), (3) DM control (diabetic, no intervention), (4) DM Gp1 (diabetic, 3 µL GH), (5) DM Gp2 (diabetic, 6 µ L GH), (6) DM Acarbose (diabetic, treated with acarbose). Fasting blood glucose levels for non-diabetic (non-DM) and diabetic (DM) groups were evaluated before and after the 10 days of diabetic induction. DM groups (excess of food and two injections of STZ) have caused a significant increment in the fasting blood glucose to 11.55 mmol/l (p < 0.0001). Both GH treatments effectively decreased postprandial blood glucose levels and the area under the curve in the oral glucose tolerance test (OSTT). Based on these results, it is concluded that green honey could play a role in hyperglycemia management and show potential as a natural alternative to conventional diabetes therapy. The underlying mechanisms need to be clarified, and their potential use in human diabetes therapy needs to be investigated.

Molecular docking and dynamics simulation studies uncover the host-pathogen protein-protein interactions in Penaeus vannamei and Vibrio parahaemolyticus.

Shrimp aquaculture contributes significantly to global economic growth, and the whiteleg shrimp, Penaeus vannamei, is a leading species in this industry. However, Vibrio parahaemolyticus infection poses a major challenge in ensuring the success of P. vannamei aquaculture. Despite its significance in this industry, the biological knowledge of its pathogenesis remains unclear. Hence, this study was conducted to identify the interaction sites and binding affinity between several immune-related proteins of P. vannamei with V. parahaemolyticus proteins associated with virulence factors. Potential interaction sites and the binding affinity between host and pathogen proteins were identified using molecular docking and dynamics (MD) simulation. The P. vannamei-V. parahaemolyticus protein-protein interaction of Complex 1 (Ferritin-HrpE/YscL family type III secretion apparatus protein), Complex 2 (Protein kinase domain-containing protein-Chemotaxis CheY protein), and Complex 3 (GPCR-Chemotaxis CheY protein) was found to interact with -4319.76, -5271.39, and -4725.57 of the docked score and the formation of intermolecular bonds at several interacting residues. The docked scores of Complex 1, Complex 2, and Complex 3 were validated using MD simulation analysis, which revealed these complexes greatly contribute to the interactions between P. vannamei and V. parahaemolyticus proteins, with binding free energies of -22.50 kJ/mol, -30.20 kJ/mol, and -26.27 kJ/mol, respectively. This finding illustrates the capability of computational approaches to search for molecular binding sites between host and pathogen, which could increase the knowledge of Vibrio spp. infection on shrimps, which then can be used to assist in the development of effective treatment.

The first ITS2 sequence data set of eDNA from honey of Malaysian giant honeybees (Apis dorsata) and stingless bees (Heterotrigona itama) reveals plant species diversity.

OBJECTIVES: Pollen is a useful tool for identifying the provenance and complex ecosystems surrounding honey production in Malaysian forests. As native key pollinators in Malaysia, Apis dorsata and Heterotrigona itama forage on various plant/pollen species to collect honey. This study aims to generate a dataset that uncovers the presence of these plant/pollen species and their relative abundance in the honey of A. dorsata and H. itama. The information gathered from this study can be used to determine the geographical and botanical origin and authenticity of the honey produced by these two species. RESULTS: Sequence data were obtained for both A. dorsata and H. itama. The raw sequence data for A. dorsata was 5 Mb, which was assembled into 5 contigs with a size of 6,098,728 bp, an N50 of 15,534, and a GC average of 57.42. Similarly, the raw sequence data for H. itama was 6.3 Mb, which was assembled into 11 contigs with a size of 7,642,048 bp, an N50 of 17,180, and a GC average of 55.38. In the honey sample of A. dorsata, we identified five different plant/pollen species, with only one of the five species exhibiting a relative abundance of less than 1%. For H. itama, we identified seven different plant/pollen species, with only three of the species exhibiting a relative abundance of less than 1%. All of the identified plant species were native to Peninsular Malaysia, especially the East Coast area of Terengganu. DATA DESCRIPTION: Our data offers valuable insights into honey's geographical and botanical origin and authenticity. Metagenomic studies could help identify the plant species that honeybees forage and provide preliminary data for researchers studying the biological development of A. dorsata and H. itama. The identification of various flowers from the eDNA of honey that are known for their medicinal properties could aid in regional honey with accurate product origin labeling, which is crucial for guaranteeing product authenticity to consumers.

Molecular interactions of trichoderma ß-1,4-glucosidase (ThBglT12) with mycelial cell wall components of phytopathogenic Macrophomina phaseolina.

Efficacy of a ß-1,4-glucosidase from Trichoderma harzianum T12 (ThBglT12) in disrupting the cell wall of the phytopathogenic fungus M. phaseolina (Macrophomina phaseolina) was studied, as the underlying molecular mechanisms of cell wall recognition remains elusive. In this study, the binding location identified by a consensus of residues predicted by COACH tool, blind docking, and multiple sequence alignment revealed that molecular recognition by ThBglT12 occurred through interactions between the α-1,3-glucan, ß-1,3-glucan, ß-1,3/1,4-glucan, and chitin components of M. phaseolina, with corresponding binding energies of -7.4, -7.6, -7.5 and -7.8 kcal/mol. The residue consensus verified the participation of Glu172, Tyr304, Trp345, Glu373, Glu430, and Trp431 in the active site pocket of ThBglT12 to bind the ligands, of which Trp345 was the common interacting residue. Root mean square deviation (RMSD), root mean square fluctuation (RMSF), total energy, and minimum distance calculation from molecular dynamics (MD) simulation further confirmed the stability and the closeness of the binding ligands into the ThBglT12 active site pocket. The h-bond occupancy by Glu373 and Trp431 instated the role of the nucleophile for substrate recognition and specificity, crucial for cleaving the ß-1,4 linkage. Further investigation showed that the proximity of Glu373 to the anomeric carbon of ß-1,3/1,4-glucan (3.5 Å) and chitin (5.5 Å) indicates the nucleophiles' readiness to form enzyme-substrate intermediates. Plus, the neighboring water molecule appeared to be correctly positioned and oriented towards the anomeric carbon to hydrolyze the ß-1,3/1,4-glucan and chitin, in less than 4.0 Å. In a nutshell, the study verified that the ThBglT12 is a good alternative fungicide to inhibit the growth of M. phaseolina.Communicated by Ramaswamy H. Sarma.

In silico assessment of dehalogenase from Bacillus thuringiensis H2 in relation to its salinity-stability and pollutants degradation.

Increased scientific interest has led to the rise in biotechnological uses of halophilic and halotolerant microbes for hypersaline wastewater bioremediation. Hence, this study performed molecular docking, molecular dynamic (MD) simulations, and validation by Molecular Mechanic Poisson-Boltzmann Surface Area (MM-PBSA) calculations on the DehH2 from Bacillus thuringiensis H2. We aimed to identify the interactions of DehH2 with substrates haloacids, haloacetates, and chlorpyrifos under extreme salinity (35% NaCl). MD simulations revealed that DehH2 preferentially degraded haloacids and haloacetates (-6.3 to -4.7 kcal/mol) by forming three or four hydrogen bonds to the catalytic triad, Asp125, Arg201, and Lys202. Conversely, chlorpyrifos was the least preferred substrate in both MD simulations and MM-PBSA calculations. MD simulation results ranked the DehH2-L-2CP complex (RMSD □0.125-0.23 nm) as the most stable while the least was the DehH2-chlorpyrifos complex (RMSD 0.32 nm; RMSF 0.0 - 0.29). The order of stability was as follows: DehH2-L-2CP > DehH2-MCA > DehH2-D-2CP > DehH2-3CP > DehH2-2,2-DCP > DehH2-2,3-DCP > DehH2-TCA > DehH2-chlorpyrifos. The MM-PBSA calculations further affirmed the DehH2-L-2CP complex's highest stability with the lowest binding energy of -45.14 kcal/mol, followed closely by DehH2-MCA (-41.21 kcal/mol), DehH2-D-2CP (-31.59 kcal/mol), DehH2-3CP (-30.75 kcal/mol), DehH2-2,2- DCP (-29.72 kcal/mol), DehH2-2,3-DCP (-22.20 kcal/mol) and DehH2-TCA (-18.46 kcal/mol). The positive binding energy of the DehH2-chlorpyrifos complex (+180.57 kcal/mol) proved the enzyme's non-preference for the substrate. The results ultimately illustrated the unique specificity of the DehH2 to degrade the above-said pollutants under a hypersaline condition.Communicated by Ramaswamy H. Sarma.

Molecular docking and molecular dynamics simulations of a mutant Acinetobacter haemolyticus alkaline-stable lipase against tributyrin.

We previously reported on a mutant lipase KV1 (Mut-LipKV1) from Acinetobacter haemolyticus which optimal pH was raised from 8.0 to 11.0 after triple substitutions of surface aspartic acid (Asp) with lysine (Lys). Herein, this study further examined the Mut-LipKV1 by molecular docking, molecular dynamics (MD) simulations and molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) calculations to explore the structural requirements that participated in the effective binding of tributyrin and its catalytic triad (Ser165, Asp259 and His289) and identify detailed changes that occurred post mutation. Mut-LipKV1 bound favorably with tributyrin (-4.1 kcal/mol) and formed a single hydrogen bond with His289, at pH 9.0. Despite the incongruent docking analysis data, results of MD simulations showed configurations of both the tributyrin-Mut-LipKV1 (RMSD 0.3 nm; RMSF 0.05 - 0.3 nm) and the tributyrin-wildtype lipase KV1 (tributyrin-LipKV1) complexes (RMSD 0.35 nm; RMSF 0.05 - 0.4 nm) being comparably stable at pH 8.0. MM-PBSA analysis indicated that van der Waals interactions made the most contribution during the molecular binding process, with the Mut-LipKV1-tributyrin complex (-44.04 kcal/mol) showing relatively lower binding energy than LipKV1-tributyrin (-43.83 kcal/mol), at pH 12.0. All tributyrin-Mut-LipKV1 complexes displayed improved binding free energies over a broader pH range from 8.0 - 12.0, as compared to LipKV1-tributyrin. Future empirical works are thus, important to validate the improved alkaline-stability of Mut-LipKV1. In a nutshell, our research offered a considerable insight for further improving the alkaline tolerance of lipases.Communicated by Ramaswamy H. Sarma.

The assessment of three dimensional modelling design for single strand DNA aptamers for computational chemistry application.

Aptamers are oligonucleotides and peptides around 15-100 bases in length and are suitable as detection probes or as therapeutics molecules. There are growing interests in the aptamer screening approach through computational simulation methods. DNA and RNA modelling lacks of validation on their predicted 3D structures due to less number of validation tools, unlike protein structures. We suggest an approach to design the stem-loop/hairpin for the three dimensional structure of DNA aptamers through serial applications of computational prediction methods by comparing the simulated structures with the experimental data deposited in PDB Data bank, followed by MD simulations. The result shows minimal structural differences were observed between the designed and the original NMR aptamers, and the stem-loop conformational structures were also retained during the MD thus suggesting the proposed aptamers designing methods are able to synthesize a high quality molecular structure of hairpin aptamers, comparable to the NMR structures.

Comparative modelling studies of fruit bromelain using molecular dynamics simulation.

Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.

In silico mutation on a mutant lipase from Acinetobacter haemolyticus towards enhancing alkaline stability.

Alkaline-stable lipases are highly valuable biocatalysts that catalyze reactions under highly basic conditions. Herein, computational predictions of lipase from Acinetobacter haemolyticus and its mutant, Mut-LipKV1 was performed to identify functionally relevant mutations that enhance pH performance under increasing basicity. Mut-LipKV1 was constructed by in silico site directed mutagenesis of several outer loop acidic residues, aspartic acid (Asp) into basic ones, lysine (Lys) at positions 51, 122 and 247, followed by simulation under extreme pH conditions (pH 8.0-pH 12.0). The energy minimized Mut-LipKV1 model exhibited good quality as shown by PROCHECK, ERRAT and Verify3D data that corresponded to 79.2, 88.82 and 89.42% in comparison to 75.2, 86.15, and 95.19% in the wild-type. Electrostatic surface potentials and charge distributions of the Mut-LipKV1 model was more stable and better adapted to conditions of elevated pHs (pH 8.0 - 10.0). Mut-LipKV1 exhibited a mixture of neutral and positive surface charge distribution compared to the predominantly negative charge in the wild-type lipase at pH 8.0. Data of molecular dynamics simulations also supported the increased alkaline-stability of Mut-LipKV1, wherein the lipase was more stable at a higher pH 9.0 (RMSD = ∼0.3 nm, RMSF = ∼0.05-0.2 nm), over the optimal pH 8.0 of the wild-type lipase (RMSD = 0.3 nm, RMSF = 0.05-0.20 nm). Thus, the adaptive strategy of replacing surface aspartic acid to lysine in lipase was successful in yielding a more alkaline-stable Mut-LipKV1 under elevated basic conditions.Communicated by Ramaswamy H. Sarma.

Substrate docking and molecular dynamic simulation for prediction of fungal enzymes from Trichoderma species-assisted extraction of nanocellulose from oil palm leaves.

Fungi of the Trichoderma species are valued industrial enzymes in support of the 'zero-waste' technology to convert agro-industrial biomass into valuable products, i.e. nanocellulose (NC). In this study, an in silico approach using substrate docking and molecular dynamic (MD) simulation was used to predict the order of which the multilayers of cellulosic polymers, i.e. lignin, hemicellulose and cellulose in oil palm leaves (OPL) are degraded by fungal enzymes, endocellulase and exocellulase. The study aimed to establish the catalytic tendencies of the enzymes to optimally degrade the cellulosic components of OPL for high yield production of NC. Energy minimized endocellulase and exocellulase models revealed satisfactory scores of PROCHECK (90.0% and 91.2%), Verify3D (97.23% and 98.85%) and ERRAT (95.24% and 91.00%) assessments. Active site prediction by blind docking, COACH meta-server and multiple sequence alignment indicated the catalytic triads for endocellulase and exocellulase were Ser116-His205-Glu249 and Ser382-Arg124-Asp385, respectively. Binding energy of endocellulase docked with hemicellulose (-6.0 kcal mol-1) was the most favourable followed by lignin (-5.6 kcal mol-1) and cellulose (-4.4 kcal mol-1). Exocellulase, contrarily, bonded favorably with lignin (-8.7 kcal mol-1), closely followed by cellulose (-8.5 kcal mol-1) and hemicellulose (-8.4 kcal mol-1). MDs simulations showed that interactions of complexes, endocellulase-hemicellulose and the exocellulase-cellulose being the most stable. Thus, the findings of the study successfully identified the specific actions of sugar-acting enzymes for NC production. Communicated by Ramaswamy H. Sarma.

Comparative structural analysis of fruit and stem bromelain from Ananas comosus.

Cysteine proteases in pineapple (Ananas comosus) plants are phytotherapeutical agents that demonstrate anti-edematous, anti-inflammatory, anti-thrombotic and fibrinolytic activities. Bromelain has been identified as an active component and as a major protease of A. comosus. Bromelain has gained wide acceptance and compliance as a phytotherapeutical drug. The proteolytic fraction of pineapple stem is termed stem bromelain, while the one presents in the fruit is known as fruit bromelain. The amino acid sequence and domain analysis of the fruit and stem bromelains demonstrated several differences and similarities of these cysteine protease family members. In addition, analysis of the modelled fruit (BAA21848) and stem (CAA08861) bromelains revealed the presence of unique properties of the predicted structures. Sequence analysis and structural prediction of stem and fruit bromelains of A. comosus along with the comparison of both structures provides a new insight on their distinct properties for industrial application.

An S188V mutation alters substrate specificity of non-stereospecific α-haloalkanoic acid dehalogenase E (DehE).

The non-stereospecific α-haloalkanoic acid dehalogenase E (DehE) degrades many halogenated compounds but is ineffective against ß-halogenated compounds such as 3-chloropropionic acid (3CP). Using molecular dynamics (MD) simulations and site-directed mutagenesis we show here that introducing the mutation S188V into DehE improves substrate specificity towards 3CP. MD simulations showed that residues W34, F37, and S188 of DehE were crucial for substrate binding. DehE showed strong binding ability for D-2-chloropropionic acid (D-2CP) and L-2-chloropropionic acid (L-2CP) but less affinity for 3CP. This reduced affinity was attributed to weak hydrogen bonding between 3CP and residue S188, as the carboxylate of 3CP forms rapidly interconverting hydrogen bonds with the backbone amide and side chain hydroxyl group of S188. By replacing S188 with a valine residue, we reduced the inter-molecular distance and stabilised bonding of the carboxylate of 3CP to hydrogens of the substrate-binding residues. Therefore, the S188V can act on 3CP, although its affinity is less strong than for D-2CP and L-2CP as assessed by Km. This successful alteration of DehE substrate specificity may promote the application of protein engineering strategies to other dehalogenases, thereby generating valuable tools for future bioremediation technologies.

Identification of functional residues essential for dehalogenation by the non-stereospecific α-haloalkanoic acid dehalogenase from Rhizobium sp. RC1.

The non-stereospecific α-haloalkanoic acid dehalogenase DehE from Rhizobium sp. RC1 catalyzes the removal of the halide from α-haloalkanoic acid D,L-stereoisomers and, by doing so, converts them into hydroxyalkanoic acid L,D-stereoisomers, respectively. DehE has been extensively studied to determine its potential to act as a bioremediation agent, but its structure/function relationship has not been characterized. For this study, we explored the functional relevance of several putative active-site amino acids by site-specific mutagenesis. Ten active-site residues were mutated individually, and the dehalogenase activity of each of the 10 resulting mutants in soluble cell lysates against D- and L-2-chloropropionic acid was assessed. Interestingly, the mutants W34→A,F37→A, and S188→A had diminished activity, suggesting that these residues are functionally relevant. Notably, the D189→N mutant had no activity, which strongly implies that it is a catalytically important residue. Given our data, we propose a dehalogenation mechanism for DehE, which is the same as that suggested for other non-stereospecific α-haloalkanoic acid dehalogenases. To the best of our knowledge, this is the first report detailing a functional aspect for DehE, and our results could help pave the way for the bioengineering of haloalkanoic acid dehalogenases with improved catalytic properties.

Insights into the stereospecificity of the d-specific dehalogenase from Rhizobium sp. RC1 toward d- and l-2-chloropropionate.

Halogenated compounds are recalcitrant environmental pollutants prevalent in agricultural fields, waste waters and industrial by-products, but they can be degraded by dehalogenase-containing microbes. Notably, 2-haloalkanoic acid dehalogenases are employed to resolve optically active chloropropionates, as exemplified by the d-specific dehalogenase from Rhizobium sp. RCI (DehD), which acts on d-2-chloropropionate but not on its l-enantiomer. The catalytic residues of this dehalogenase responsible for its affinity toward d-2-chloropropionate have not been experimentally determined, although its three-dimensional crystal structure has been solved. For this study, we performed in silico docking and molecular dynamic simulations of complexes formed by this dehalogenase and d- or l-2-chloropropionate. Arg134 of the enzyme plays the key role in the stereospecific binding and Arg16 is in a position that would allow it to activate a water molecule for hydrolytic attack on the d-2-chloropropionate chiral carbon for release of the halide ion to yield l-2-hydroxypropionate. We propose that within the DehD active site, the NH group of Arg134 can form a hydrogen bond with the carboxylate of d-2-chloropropionate with a strength of ∼4 kcal/mol that may act as an acid-base catalyst, whereas, when l-2-chloropropionate is present, this bond cannot be formed. The significance of the present work is vital for rational design of this dehalogenase in order to confirm the involvement of Arg16 and Arg134 residues implicated in hydrolysis and binding of d-2-chloropropionate in the active site of d-specific dehalogenase from Rhizobium sp. RC1.